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Abstract. Consider a droplet of liquid on top of a grooved substrate. The wetting or 
not of a groove implies the crossing of a potential barrier as the interface has to distort, 
to hit the bottom of the groove. We start with computing the free energies of the dry 
and wet states in the context of a simple thermodynamical model before switching to 
a random microscopic version pertaining to the Solid-on-Solid (SOS) model. For some 
range in parameter space (Young angle, pressure difference, aspect ratio), the dry and 
wet states both share the same free energy, which means coexistence. We compute these 
coexistence lines together with the metastable regions. In the SOS case, we describe the 
dynamic transition between coexisting states in wetting. We show that the expected 
time to switch from one state to the other grows exponentially with the free energy bar- 
rier between the stable states and the saddle state, proportional to the groove's width. 
This random time appears to have an exponential-like distribution. 

PACS classification: 68.08.Bc, 68.03.Cd , 05.40.-a. 
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1. Introduction and outline 

The detailed study of superhydrophobic surfaces has revealed that on a rough substrate, 
a drop can present two shapes: either one obeying the Cassie-Baxter equation or one obey- 
ing the Wenzel equation [6] , [10] , [5] . On the other hand, any transition between two states 
depends on the height of the barrier which has to be overcome. The corresponding transi- 
tion will thus be a function of time as revealed already by Kramers law. We review here 
this crucial aspect of the problem within the framework of an exactly solvable statistical 
mechanics model. 

The corresponding Cassie-Baxter and Wenzel states are illustrated in Fig. [T] In the 
first (dry) state, there is no wetting of the sides of the U-shaped wells, with vapor trapped 
in-between, whereas in the second (wet) one, there is at least a partial wetting of the bottom 
of the U, with vapor trapped in the corners. This is illustrated in Fig. [21 zooming on a 
single well. The transition between these two states is related to the amplitude of some free- 
energy barrier: in this rough picture, this barrier can be overcome by applying an external 
force triggered for instance by an impact velocity or equivalently by increasing the pressure 
difference Ap between liquid and vapor. A classical way to characterize the wet shape is 
given by the Young contact angle 9, as sketched on Fig. (3] 



A rough substrate is thus made of clefts (wells, grooves...) governing wetting properties. 
For a better understanding of this physical problem, we need first to analyze and characterize 
the free energies of both states as functions of the parameters 0, Ap and p, the latter being 
the shape factor of a typical well. We first do that in Section 2 for two simple macroscopic 
solvable models of wetting, one simpler without the pressure parameter and the other one 
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including pressure. For the isotropic model including pressure, we show that there exists a 
region in parameter space where the two states can both exist (the free energy function has 
two stable minima) and therefore regions where only a single of these two states exists (the 
free energy function possesses a single stable minimum). Within the metastable region of 
existence of the two states, we exhibit a line of coexistence along which the free energies of 
both states are exactly equal. 

We can easily imagine that, due to fluctuations, there will be an opportunity to flip be- 
tween the 'dry' and 'wet' states of Fig. [21 So we need to embed our wetting problem into 
the framework of a stochastic model of wetting. For this purpose, we used the microscopic 
statistical mechanics SOS model. In this context also, the wetting of a well implies the 
crossing of a potential barrier: the interface has to distort, to hit the bottom of the well. 
We wish to understand this problem in more details. 

The aim of Section 3 is to show that we are still able to compute analytically the equi- 
librium free energies of the dry and wet phases in this statistical SOS context. Let us be 
more precise. 

We study our wetting problem in the context of the SOS model in a square well of width 
I —: nlo and height /i, with h/l ~: p. The equilibrium measure for the interface heights 
/li > is given by: 

n 

(1.1) d/.„ = ^-ie-('^S;UI'«.+i-'«.|+f EIU/^O J|(l + a5(/i,))d/i„ 

with /iQ = hn+i ~ h. Here Z„ is the normalizing partition function, and the pressure 
term comes from Ap h{x)dx ~ ^P^i^i-^ = ^J^^i^ which defines K a.s K := lAp. 
The length Iq = l/n is the coarse grain length at which the SOS model is defined as an 
approximation to a truly microscopic model. Thanks to this coarse-graining, the interface, 
originally of finite- width, constrained by the corners of the well, has become a surface pinned 
at the corners. 



The SOS model has three positive parameters (J, K, a), representative for the first two of 
the length of the interface and the area below the interface, and for the third, of the liquid's 
affinity for the bottom of the well. More precisely, 

dlogZn /^., ,,\ , dlogZ„ ,/^,\ , d\ogZ„ /" \ 
dJ- ^ 1^'+^ - log^^ = lo {^h.^ , /o^^ ^ ° I 

are respectively the mean excess length of interface, the mean area below the interface, and 
the mean wetted length of well. These three parameters are in correspondence with the 
surface tension aL,v between the liquid and vapor phases, the pressure difference Ap and 
the Young angle 0. 

In the thermodynamic limit, the structure of the interfaces is given by a Wulff shape 
pinned at the corners of the well, expressed analytically in terms of the projected surface 
tension. Two cases arise (see Fig. [2]): 

- 'dry' case: the Wulff shape does not hit the bottom of the well and so hangs between 
the two corners. 

- 'wet' case: the Wulff shape does hit the bottom of the well and the interface is made of 
three pieces whose central part, fiat (corresponding to the wetting of the substrate by the 
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Figure 1 . The Cassie- Wenzel states 



cassiewenz 




Dry 



Wet 



V 



Figure 2. Dry versus Wet states 
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fluid), is linked up to the corners of the well by two symmetric pieces of the Wulff type. 

For the sake of simplicity, in this study, the vertical sides of the well are completely 
hydrophobic. This fits with the SOS model (|l.ip with hi eM. for which the Wulff shape has 
no vertical part. 

Using the results obtained in [l]-[4], [8]- [9], we can compute exactly the normalized free 
energy for each of these two configurations. At fixed J and p, we compute the line of coexis- 
tence a = a (K) for which these free energies coincide. Due to a one-to-one correspondence 
between a and 9, a phase diagram 9 = 6 (K) follows. In this phase diagram, we compute 
the lines separating a metastable region where the two phases coexist and stable regions 
where a single phase (either wet or dry) is stable. 

When both 'wet' and 'dry' SOS states share the same free energies, we are subsequently 
interested by the transition between them. In the Monte-Carlo dynamics of Section 4, we 
show that the system undergoes rare transitions between these two equilibrium states as time 
passes by. We also study the first passage times from one state to the other, together with 
finite-size corrections. This requires some preliminary understanding of the free energy of an 
unstable SOS saddle state which can be computed explicitly, in the thermodynamic limit. 
We show that the expected time to switch from one state to the other grows exponentially 
with the free energy barrier between the stable states and the saddle state, proportional to 
the system's size. This random time turns out to have an exponential-like distribution. 

In Appendix C, a toy Markov-chain model is designed which illustrates the typical be- 
haviors encountered in our random wetting problem. The rate of growth of the expected 
time to move from one stable state to the other when an unstable state lies in-between is 
the energy barrier between the unstable and stable states (Kramers' law), and the random 



4 



JOEL DE CONINCK FRANCOIS DUNLOP THIERRY HUILLET 



time itself normalized by its mean converges in distribution to an exponential probability 
distribution. 

2. A SIMPLE MACROSCOPIC MODEL. EQUILIBRIUM FREE ENERGIES 

Before we run into similar considerations pertaining to the SOS model, let us first briefly 
describe some easy macroscopic arguments that led us to the forthcoming study. 

2.1. No pressure {K = 0) 

. DRY PHASE. Consider a well of length / and height h = pi. Call it the substrate (S). Fix 
the origin of the {x, z) axis at the middle of the bottom of the well. In the following, we 
take / as unit of length, so that the bottom of the well is the interval [—0.5, +0.5] at height 
2 = 0. 

The well is initially filled with gas or vapor {V). We wish to fill this groove with a liquid 
(L) in the absence of pressure (Ap = f = 0). We assume that the vertical walls of the well 
remain dry. In a dry phase, there exists a Wulff shape separating the LV phases which is 
pinned at {x = —0.5, z = p) and {x = 0.5, z = p) without hitting the bottom of the well: it 
is the straight line joining these two corner points. The well is entirely filled with vapor, 
the liquid entirely stands above this straight separating line and so the liquid does not wet 
the substrate at all. The specific free energy of this dry phase is 

S (2.1) /drv := ^ = ^Ly + (1 + 2p) • ays, 

where gab is the surface tension between the phases A and B. 

WET PHASE [K = 0). In a wet phase, the liquid will meet the substrate at the bottom 
of the well. Let 9 be the Young angle between the vapor phase and the substrate at the 
right-most meeting point. The larger 6 is, the more the substrate is hydrophilic. 

Let xi E (0,0.5). In a wet phase with K = 0, the WulfF shape would indeed be made 
of two symmetric linear pieces of length separating phase L from V, one joining 

point (x = —0.5, z = p) to point {x = —Xi, z = 0) and the other joining {x = x\,z = 0) to 
(x = 0.5, z = p) . In between, the line joining (x = —xi, z = 0) to {x = xi, z = 0) would be 
the flat wetting zone for the substrate. For a given 6, a wet phase can exist if and only if 
a;i > or else 

"12] (2.2) p< 0.5 tan 6* 

(the well is not too deep). 

When this wet phase exists, its specific free energy is 

Is] (2.3) Uet~^=2^-aLv + {l + 2p-2xi)-avs + '2xi-aLS- 

' ' I smO 

When both phases exist, the question of which phases is favored makes sense. So we can 

ask for conditions under which F^ry < F^et- Using aLS — (^vs = CLycos^, we easily get 

the condition 

"14] (2.4) p> 0.5tan(6»/2), 

meaning that for the dry phase to win over the wet one, the depth of the well has to be 
large enough. 

So, when a wet phase exists (p < 0.5 tan 0), the dry phase wins over the wet phase 
whenever the well's depth satisfies p > 0.5 tan (0/2). Else, if p > 0.5 tan 6, this problem does 
not make sense simply because only the dry phase exists and so it necessarily wins. So, in 
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this oversimplified wetting problem, we expect a parameter range for which the two phases 
wet and dry can coexist (a metastable region), and within this parameter range another 
parameter range for which one phase is favorable over the other. When p = 0.5 tan (6'/2), 
the two phases share the same free energy (a bistable line of coexistence of both phases). 

2.2. Pressure [K > 0). WET PHASE. The Wulff-shaped lines separating the LV phases 
are now symmetric convex circle arcs with radius 

''^ Ap ~ K - 

Define the scaled radius as R — r/l. In the scaled length unit and in the wet phase, we 
have one circle arc joining point (x = —0.5, z = p) to point {x = —Xi, z = 0) and the other 
joining A := {x = Xi, z — 0) to B :— {x — 0.5, z ~ p) . 

Let O be the center of the latter circle. Let be the angle {OA, OB) . Then the scaled 
euclidean distance between A and B is given by (see Appendix Al) 

AB^Ur' (sinH+P^-VA''^' 



R 



where A = (sin^ 9 + - {-^f. Next AB = 2R avcsin {AB / {2R)) is the scaled arc 

length of the arc joining A to B. Using 

AB = ^J{0.5- xi f + p^. 
gives xi as an explicit function of {0, p, R). 

For a given 9, the wet phase can exist if and only if xi > which is, observing tan {9 + Lp/2) = 
p 



0.5-a:i ■ 

p < 0.5 tan {9 + (f/ 2) 
(again, the well should not be too deep). 

When this wet phase exists, its specific free energy is found to be 
lb] (2.5) U,t-=^ = '2AB-aLv + {l + 2p-2xi)-avs + 2xi-aLs + 2K-A^,t, 

where A^jj^t = f (0.5 — xi) — i [R"^^ — R^ smip] is the scaled dry area of the vapor beneath 
the right part of the Wulff shape. The right-most pressure term in (|2.5p comes from 

j{Ap-l^A^et) =K-A^et. 

DRY PHASE. With bIb = 2i?arcsin {j^) =: R^), the arc length between the left and 
right corner points B' and B of the well, f^et should be compared to the specific free energy 
in the dry phase which is 

ji] (2.6) fdry.^^=B^-^LV + {l + '2p)-CTvS+KAdry. 

Here Adry ~ p—\ [R^i^ — R^ sin -0] is the scaled dry area of the vapor beneath the hanging 
Wulff line anchored at B' and B. 

This leads again to an implicit critical line of coexistence in the parameter space where 
fdry — fwet- In Fig. [31 wc plot the line of coexistence 9 = 9 (K) for the arbitrary set 
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Figure 3. The line of coexistence F^gt — Fdry from (|2.5|l2.6p with aLV — 
(TVS = 0.01 and p = 0.2. 



of parameters aj^y — ays — 0.01 and p = 0.2. The largest value of 6 is obtained when 
if = as 6'inax — 2 arctan (p/0.5) , from (|2.4I) . The largest value of K is obtained from the 
dry phase when the circle arc pinned at the corners of the U well hits the bottom of the 
substrate tangentially (9 — 0). We get i^Tmax = 2paLv/ {p^ + 0.5^) . 

3. SOS MODEL: EQUILIBRIUM FREE ENERGIES 

We now run into similar considerations for the SOS model of wetting arising in statistical 
mechanics. In this Section, we compute the free energies of the dry and wet phases in the 
thermodynamic limit for the SOS model ()1.1|) . For the sake of simplicity, we decided to 
work at fixed values of the parameters J — 3.0 and p = 0.2, leaving the free parameter 
space be restricted to a and K. 

DRY PHASE. In the context of the SOS model (fTT]) . the typical Wulff shapes which 
come in can be described as follows. Let 



with derivative 



fit) = {i + jH') 
fit)- 



1/2 



1 



/(t) + l 

The projected interface tension aLv corresponding to (jl.ll) with /i^ g R and a = and 
boundary conditions implying a slope tan^ = t, with energy measured in units of kT, is 
defined as 



(3.1) 



CTLV (t) = \im log Z„ {t) a (t) , 

n— >oo n 
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where 

'fit) 



M (3-2) 5(t) = /(t)-log, J 

We have a' (t) = f([-^^2 • Then the relevant Wulff shape equations relative to (ll.ip and (|3.1 
are implicitly given by (see [4], |8] and [9]) 

Z it) ^ --lo i^Lv (t) ~ ta'^v (t)) ; Z (0) = -l„ log i^-j 
where t is the tangent dZ/dX to the curve X ^ Z (X) . 

Introduce the scaled variables x (t) = X (t) /I and z (t) = Z (t) /I. Then, recalling I — nl 

"wTl (3.3) x{t) = (t); x{0) ^0 
K 

[W2] (3.4) z (t) = -1 (5 it) - ta' {t)) ; H0) = ^ log 

where the range of x (t) is (— ■^) and where t is the tangent dz/dx to the curve x — > z (x) , 







Note that t (x) 



2Kx 



The scaled well has now unit length and fixed height p = 0.2. As before, fix the origin of 
the [x, z) axis at the middle of the bottom of the well. We wish to derive the equation of 
a Wulff shape p. 31 13. 4p which is pinned at {x = —0.5, z = p) and {x = 0.5, z = p) without 
hitting the bottom of the well. Let to be the value of the tangent at {x = 0.5, z — p), with 
to = t (0.5) = pJk'^ /4 - have 

Zdry (t) = z{t) - Z (to) + p 

SO that Zdry (0) ^ z{0)- z (to) + 0.2 = log (f ) - z (to) + p > 0. 

When K increases, the minimum of the hanging Wulff shape gets closer to the bottom of 
the well. There is a value K = -flTmax for which this minimum hits tangentially the bottom 
of the well in one point. The value of if max is characterized by 

.... (0) k.K.„..= ^ log (§) - z {jT^^) +P-0. 

The admissible range of K in the dry regime is thus [0, i^max] ■ 

In Appendix A2a, we obtain the specific free energy of the dry phase as 

' 1 + 0.5/b " 
1 - 0.5/6 ' 

Note that when K 0, to and z (to) ~ z (0) = log ) . Thus fdry log (j) = a (0) 
which is the specific free energy of the trivial {K = 0) —Wulff shape joining linearly the cor- 
ners of the well parallel to the bottom of the well. 

WET PHASE. Let xi e (0,0.5). In the wet phase, the Wulff shape is made of two 
symmetric convex pieces, one joining point {x — —0.5, z = p) to point {x = —xi, z = 0) and 



1^ (3.5) fdry = K{0.2~z{to))-2 + 2blog 
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the other jommg (x — xi, z — 0) to {x — 0.5, z — p) . In between, from (x = —xi, 2: = 0) to 
(a; = xi, z = 0), the curve is flat, pinned to the substrate. 
From [5] , the specific free energy of the wet part is 

swet I (3.6) o-u,et = -log I -^^-^ ) , a > flc := 1/J. 

' ' \ZaJ — 1/ 

Consider the right part of the Wulff shape. Let ti be the tangent of the Young angle 9, 
which is also the interface slope at {x = xi, z — 0); assuming a > a^, we have (see (5) (8) (9) 

mm 

I 1 /„ _N o 2a(aJ—l) 

Young (3.7) h = tane* = — ^ > 0. 

I 1 2a J — 1 

Note that when a < Qc, the affinity for the bottom of the well is not strong enough to 
produce a wet part in the equilibrium Wulff shape. 

Let t2 be the tangent of the angle of the Wulff shape at point {x = 0.5, z — p) . Using the 
canonical equation p. 3) 13. 4|) of a standard Wulff shape, the equations of this Wulff shape 
are given by: 

Xwet (t) ^ x{t) - X (^2) + 0.5 
Zwet {t) = Z (t) - Z (ti) , 

where t2 has to be determined implicitly by z^et (^2) = P- We then have 

xi = Xwet {h) = X {ti) - X (12) + 0.5. 

The specific free energy of the wet phase is obtained as (see Appendix A2b) 



Fwet (3.8) /^et = 2a;i(7„et-2(0.5-a;i)(is:z(ti) + 2) + 261og 



^-^ ^ xi+a:(t2)-0.5 ^ _ £(£2)^ 



FwetKO 



When K — 0, the two pieces of the Wulff shapes become straight lines. The tangent of 
the Young angle is thus <i = p/ (0.5 — xi) and 



(3.9) 



fwet I 



wet IR'=0- 



2 (0.5 - xi) a (ti) + 2xi(Twet- 



There exists a maximal value Cmax of a characterized by: f^et 1^=0= fdry \k=q^ log (J/2) . 
For a > Cniax, the wet phase has a lower free energy than the dry phase, for all K > 0. 

We fix J = 3.0 and look for the values a = a (K) for which = fdry, using (|3.5p and 
(|3.8p . meaning coexistence of the two phases. In this example, the range of K is [0,i4rmax] 
with Kmax — 4.6767 and the range of a is [ac,amax] with Oc — 1/J and Omax — 1.5114. 
Using this curve a — a (K), together with p.7p . relating the Young angle 6 to a, we rather 
consider the line of coexistence = [K) . This line of coexistence is shown on Fig. In 
this phase diagram plot, the dotted lines separate a metastable region where the two phases 
coexist and stable regions where a single phase (either wet or dry) is stable; the solid line of 
coexistence separates the two stable phases within the metastable region. The two dotted 
lines are obtained while using fgad — fdry and fsad = fwet, respectively. Note that the point 
dX K = Q separating the dry stable zone from the metastable zone is exactly characterized 
by p = 0.2 = 0.5 tan 61 as in 
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Figure 4. The line of coexistence fdry — fwet from <\3.5\ 13. 8|) . together 
with the metastabihty hnes: fdry = fsad from ((XSl li?T|) and /u,et = /^ad 
from dsn mil); here J = 3.0 and p = 0.2. 
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4. Dynamics and numerical simulations 

Let us first compute the equihbrium free energy of the saddle-point phase, characterized 
by a single contact point in the center of the well. 

SADDLE-POINT PHASE. In the saddle-point phase, the Wulff shape is made of two 
symmetric pieces and xi = so that there is no flat part corresponding to wetting. With 
ti the tangent of the Young angle at point (cc = 0, z = 0), the equations of the WulfF shape 
in the saddle-point configuration are 

Xsad (t) ^ x{t) - X (ii) 
Zsad (t) = z{t) - Z (ti) . 

In Appendix A2c, we show that the specific free energy of the saddle phase reads 



Fsad I (4.1) fsad = -Kz{ti)-2 + 2b\og 



which is implicitly known because so is ti 



^ ^ 0.5+a:(ti) ^ ^-y _ x{ti 



(ti) \ /-, _ 0.5+a(fi) 



1 r 



4.1. Dynamics. For a square well of width I = uIq, we consider a Markovian dynamics of 
Monte-Carlo type having (jl.ll) as invariant measure. The free energy barrier to cross starting 
from the dry (wet) phase is F^"l - f]"^ - nloifsad - fdry) (respectively F^^d " ^wlt 
nlQ{fsad — fwet))- Conventional wisdom suggests that, with {t^J2/ mean time needed to 
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first enter the wet (dry) phase starting from the dry (wet) phase, in a system of size I = nlo 

as n CXI, 



1 



:;^^0g{T''^'2} -> fsad- f dry 



1 



log {TwI) fsad - fwet- 

The limiting quantities coincide, assuming coexistence, 9 = Q{K). We expect polynomial 
corrections such as e'^'-^'^f"'^'^^^'^''^^ (r^^l^ ^ Cdwn^ and e""'°^''^=°'^"'^""") (^r^j^ ^ Cwdfi^ 
where Cdw and c^d would be two possibly distinct n— independent constants. In any case, 
the relative weights of the dry and wet phases may differ, 

' dw , \' wd 



lim (4.2) lim — r , . ^ lim 



' dw / ~^ \ ' wd / \ ' dw / ~^ \ ' wd 



4.2. Numerical simulation. Simulations were performed for one-dimensional interfaces 
over a trough of length I = n and depth h = n/5. The interface is pinned at both ends, 
ho = hn+i = n/5, and h = (hi, . . . , ft,„) e i?" is distributed at equilibrium according to 
(ini), or 

n n n 

(4.3) /x„(dh) = exp(- -h,\--Y,hi) J|(l + a5{hi))dh, 

i=0 1=1 i=l 

where the partition function Z„ normalizes the probability. 

The sub-lattice parallel heat bath dynamics, irreducible and satisfying the detailed bal- 
ance condition with respect to is defined as follows, for t = 1, 2, . . . : 

(4.4) P(dh*|h*-i) = M„(dh*) n Sihl-h^'')/ J ^,Jdh) n S{h^-fA~') 

i+t odd i+t odd 

Fig. [S] shows one hundred samples obtained from this dynamics. The interface is typically 
near one of the two Wulff shapes, "dry" or "wet" . The parameters J, K, a were chosen so 
that the free energies computed from the two Wulff shapes were approximately equal, so 
that the interface spent approximately equal times near these two shapes. The empty region 
between the two shapes indicates a region of small probability, which we shall call a free 
energy barrier. The aim of the simulation is to find the law of the escape time from dry 
to wet or conversely. The interface can be within typical fluctuation from one of the Wulff 
shapes, in which case it is easy to decide whether it is "dry" or "wet" , but large deviations 
in between cannot be attached to one or the other in any justified way. 

We measure time either in Monte-Carlo Steps per Site (MCS/S), or in Monte-Carlo 
Steps per Site divided by (MCS/S/n^), because the relaxation time of a fiat interface 
without free energy barrier is of order n? MCS/S in non-conservative dynamics. One MCS/S 
corresponds to two time steps of the form (|4.4p . This observation may well have interesting 
consequences to analyze the dynamics of spreading of nanodrops on top of rough substrates 
using molecular dynamics simulations. 

At time intervals of the order of a few MCS/S, three measurements are taken, corre- 
sponding to the observables that make up the Hamiltonian in (|4.3p : the normalized length 
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of the interface 

1 " 

(4.5) L{h) = —-Y,\h^+l-h.\, 

" i=0 

the normahzed area below the interface 

(4.6) ^(h)-^E^- 

j=i 

and the number of zeroes of h divided by 2n, denoted VF(h) (number of contacts with 
Wall). Each measurement appears as a dot on Fig. IH making up clouds of points around 
the corresponding mean values for each of the two WulfF shapes. The factors of 2 and 3 in 
the definitions of W{h) and A{h) are designed to facilitate the reading in Fig. [51 We then 
decide somewhat arbitrarily intervals of values of the three observables associated to the 
wet state or to the dry state or to the region in between, called the barrier: 

Dry : L{h) < 0.48 , A(h) > 0.27 , W{h) = 
Wet : L{h) > 0.52 , A(h) < 0.18 , W(h) > 0.04 
Barrier : otherwise 

At any given time, if all three observables agree for dry, or for wet, then the interface is 
declared dry (D), or wet (W). Otherwise it is considered to be in the barrier (B) between 
dry and wet. The barrier is clearly in a region of large deviations from the dry state or from 
the wet state. In principle one observable should be enough, but the numerical simulation 
is done with a limited set of values of n, and the picture emerges more clearly using three 
observables. 

We thus obtain a marginal of the interface dynamics, which is a non-Markovian process 
with values in {D,B,W}. Of course the simulation uses discrete time, but the mean sojourn 
time in D or B or W is of order MCS/S, so that a continuous time description is better, 
with numerical results expressed in MCS/S/n^. 
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Figure 6. From top to bottom: L{h), A{h), W{h.), function of time mea- 
sured in MCS/S/n2 (dots), with n = 700, J = 3.0, K = 4.0, a = 0.455. The 
regions between 0.48 and 0.52 for L(h), between 0.18 and 0.27 for v4(h), or 
between and 0.04 for M^(h), correspond to the barrier between 'dry' and 
'wet'. 



A marginal of this marginal is the sequence of letters, forgetting their duration, looking 
like 

DBDBDBDBWBWBWBWBWBWBWBWBDBDBDBDBDBWBWB. . . 

with only four 2-letter patterns: DB, BD, WB, BW. This restriction allows only six 
3-letter patterns: DBD, BDB, WBW, BWB, DBW, WBD. The last two are very rare 
(~ exp(— const. n)). The numerical algorithm does not strictly forbid the 2-letter patterns 
DW or WD, crossing the barrier in a few MCS/S, but they are so rare, with a relative 
frequency expected ~ exp(— const. n^), that we haven't seen them in the experiment. 

The limiting relative frequencies of the letters D, B, W, are respectively 1/4, 1/2, 1/4. 
The barrier B should be decomposed into a dry side and a wet side of the saddle point. We 
cannot describe precisely the corresponding configurations, but we can assume that within 
a 3-letter pattern DBD, the system remains in the dry state, or on the dry side of the 
saddle point, and analogously for WBW. If we would change every B in DBD into D, and 
every B in WBW into W, then we would find relative frequencies tending to 1/2, 1/2 for D 
and W. The corresponding partition of the configuration space into D and W is analogous 
to partitions which play a role in some studies of metastability [7]. There is of course a 
remainder ~ exp(— const. n) for B from DBW and WBD. 

The dynamics is started at time to = in D, with the interface near the Wulff shape 
associated with the dry state. The first time in W is denoted ti, measured in MCS/S/n^. 
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The next time in D is denoted t2, etc. Thus t2k+i is the first time in W after t2k for 
fc = 0, 1, . . . , fcmax and t2k is the first time in D after t2k-i for fc = 1, 2, . . . , fcmax- We also 
define tjfc as the last time in D before i2fe+i and after t2k, and i2fc+i ^^^^ time in 

W before t2k+2 and after t2k+i- The empirical mean escape times from D and from W are 
defined respectively as 

^ fcmax -1 ^ fcmax -1 

(4.7) TD = T J2 (*2fc - t2k) , fw = T J2 (*2fc+l " ^2fc+l) 

"-max ■ „ "'max ■ „ 

A:— A;— 

and the empirical mean barrier crossing time is defined as 

(4.8) ^B = ^ £ iti+i~t',) 

The total time of the experiment is fcmax(T^_D + 2fcmaxT^B, in MCS/S/n^. The measured 
escape times from D and from W are respectively — tj^. — t2k and = i2fc+i ~ 
t2k+i- They are not a priori independent nor identically distributed, but one may study 
the empirical cumulative distribution function of a random variable or having the 
sequence of realizations or r^: 

(4.9) R^it) = J-\{k:T^< t}\ , R'^it) = -^|{fc : < t}\ 

"-max "'max 

This is shown on Fig. [71 showing that and are approximately independent and 
identically distributed (iid) exponential random variables. The fit is better with a small 
shift as indicated, which we interpret as a redistribution of the barrier time 2fc,„ax7^s to the 
dry and wet states, in a proportion which we cannot decide directly. 

We now turn to the dependence of r^, and upon n. Fig. [8] shows that log(f-°)/n 
and log(r^)/n tend to a constant ~ 0.0105 as n — >■ oo, to be compared with the free energy 
density difference between the barrier or saddle point and the dry or wet states, giving a 
theoretical value 0.0129, cf. (|4.1l)(|3.5|)p.8p . On the other hand f^, measured like the other 
times in MCS/S/n^ appears to converge to a limit independent of n, = 0.22 ± 0.03 for 
n = 500,600,700,800. 

For each value of n, the simulations were done with a value a„ of the parameter a chosen 
so that the system spent equal time in 'dry' and in 'wet'. Finite size corrections to the 
thermodynamic limit should imply a„ = a + 0(l/n); a fit by a„ = a^o + const. /n gives 
Ooo = 4.61, with an uncertainty compatible with the expected limit a ~ 4.64, as can be read 
on Fig. Hat if = 4.0. 

The size of samples used in Fig. [8]for n = 200,300,400,500,600,700,800 was fcmax = 
15702,4936,2433,778,272,47, 19 respectively, where fcmax is the number of observed tran- 
sitions from 'dry' to 'wet', also equal to the number of observed transitions from 'wet' to 
'dry'. 



5. Conclusion 

We describe a situation where a liquid droplet is on top of a structured substrate pre- 
senting grooves or wells. We show that there exists a range of parameters for which wetting 
and non-wetting states both share the same free energy, entailing that there can be tran- 
sitions between these two states. We calculate analytically the involved free energies of 
the dry and wet states in the context first of a simple thermodynamical model and then 
for a random microscopic model based on the Solid-on-Solid (SOS) model. In the latter 
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Figure 7. Empirical cumulative distribution function of (dry) and 
(wet) and fit R°{t) = 1 - exp(-(i + 1.2)/27.8), {t) = 1 - cxp(-(t + 
1.3)/27.5), functions of time in MCS/S/n^, for n = 500, J = 3.0, = 
4.0, a = 0.451, and fc^ax = 93. 



case, we present a dynamical model showing stochastic transitions between the wetting and 
non-wetting states. We show that the expected time to switch from one state to the other 
grows exponentially with the free energy barrier between the stable states and the saddle 
state, proportional to the width of the grooves. This random time is shown to have an 
exponential-like distribution. Our study hopefully contributes to a better comprehension of 
the behavior of fluids on structured surfaces. 

6. Appendix 

Al. Proof of 

Let O be the center of the circle anchored at A and B with scaled radius R. With ip the 
angle {OA, OB) , the scaled euclidean distance between A and B is 

AB = 2i?sin (|) = \J {Q.^ - xxf + . 
The scaled arc length of the arc joining A to S is thus 

2R 



AB ^ Rip ^2R arcsin 
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for n = 200, 300, 400, 500, 600, 700, 800 respectively. 



dbw 



Using sin (6* + 1^9/2) 



AB 



and AB — 2i?sin(^) gives AB^ after eliminating ip. With A 



e + - we get AS as the exphcit fonction of {6, p, R) : 

pcos9 



AB' 



2i?2 sin^ I 



R 



Note that, geometrically, p < Rcosd which entails that A > 0. 

As a result, AB = 2i? arcsin (AS/ (2i?)) , xi = O.b-y^AB'^ ~ p"^ and ip^ 2 aicsm {AB/ (2i?)) 
have themselves an explicit expression in terms of {9,p,R). These quantities are the ones 
needed to compute fwet from (|2.5p . 



A2. Proofs of ([331), ^ and (liUl . 

A2a (DRY PHASE). The specific free energy /dry (a^i, xa) := limj^oo ^"""^f ''"""^ of the 
dry phase between {xi,X2) has two contributions, one pertaining to the length and the other 
to the area below the interface, namely: 

rX2 /•X2 /•X2 

Fdryl (6.1) fdry {xi,X2)= / adx + K I Zdrydx — / ta' dx + K {x2 — Xi) {0 .2 — Z (tf))) ■ 

J Xi J xi J xi 
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With b := J/K, we have 



ta'dx ~ K I txdx ~ 



2^2 



rdx = 2 



-2 (X2 - Xi) + 



xi 

X2 



J2 - K^X^ 

1 1 



2:2 ^2 /l2 



-2 (a;2 - xi) + 6 log 



1 + X2/6 1 - xi/b 



dx 



dx 



1 + Xi/bl- X2/b_ 

Finally, letting xi = — 0.5,X2 = 0.5, we obtain the specific free energy of the dry phase as 
in (ESI. 



A2b (WET PHASE). The limiting specific free energy of the wet phase has two parts, 
one corresponding to the symmetric pieces of the Wulff shape, the other to the flat part: 

rO.5 



(6.2) 
We have 



/•U.O 

fwet = 2 / ((? + Kz-uiet) dx^et + 2x1(7^ 



{a + Kzwet) dxwet = -Kz (ti) {X2 - Xi) + 



xi+x(t2)-0.5 



(cr + Kz) dx 



and, using 



{a + Kz)dx^ ta dx ^ -2 {l3 - a) + b\o^ 



(l + f) 

(1 + f) (l-l)' 



we finally obtain 

fwet = 2xiawet - 'i'Kz (ti) (0.5 - xi) + 2 
leading to ([3 



xit2) 



(a + Kz) dx. 



xi+x{t2)-0.5 



A2c (SADDLE PHASE). The slope t2 of the interface at point {x = 0.5, z = p) has to be 
determined implicitly by Zgad (^2) — P- We have 

Xsad (h) = X (<2) - X (ti) = 0.5 
Zsad (^2) = Z (t2) - Z {h) = p. 



From the first equation, recalling x (t) — ^ o (t) and a' (<) — 



f{t)+2 (l + ,72t2)l/2^j 



t2 = 5'-Wa' (ii) + 0.5X) =:^{ti) 



2s 



where the inverse of a' is easily seen to be 

a'-' is) = 

Thus t2 = ip (ti) is an explicit known function of ti. Plugging this expression in the second 
equation and recalling z (t) = — (a (t) — ta (i)) /K, 

{a - ta) (if (ii)) -{a- ta) (ti) = ~pK 
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giving ti implicitly and then t2 using t2 — ^ (ti). The specific free energy of the saddle-point 
phase is thus 

M.5 

Fsadl I (6.3) fsad = 2 / (ct + KZsad) dXsad- 

Jo 

We have 

{a + Kzsad) dXsad = ^Kz {ti} {X2 - xi) + j (ct + Kz) dx 

and so 

f0.5+x{ti) 



fsad = -Kz {h) + 2 {a + Kz) dx, 

Jx(ti) 

leading to (|4?T|) . 



B. A toy model. 

Although the problems encountered in this study are far from being Markovian, we find 
it useful to end up with recalling similar issues in the context of Markov chains or the like. 

Consider a discrete-time k Markov chain Xk with five states {0, 1,2,3,4}. Suppose the 
following transition probabilities Pij from state i to j hold: Po,o = 1 — a', Po,i = ct'', 
Pi,o = 1 - a, Fi,2 = a; ^2,1 = 1/2, P2,3 = 1/2; ^3,2 = a, ^3,4 = 1 - a; P^^s = a', 
P4 4 = 1 - a'. 

The parameters a and a' are small, with 

a = e"['^(2)"'^(^)l/^ = g-[c/(2)-c/(i)]/e 

the energy barrier terms within the brackets being all positive and e small. Thus {0} and 
{4} are two stable states separated by a barrier state {2}. Let us compute the law of the 
time To, 4 = inf (A; : Xk = 4 | Xq = 0) needed to move from state {0} to state {4} . The chain 
is a nearest neighbors birth and death chain which is ergodic. Putting p^, = Px,x+i and 
qx — Px,x-i, the invariant measure is ttx = ttq 0^=0 'q^- Starting from {0}, the sample 
paths are made of iid excursions separating consecutive visits to {0} . The law of the height 
H of an excursion is given by 

ip{h) 

where 

x-l y 

is the scale function of the chain. In particular, we get Pr {H = 4) = l/(p (4) = a/2. 
With fi the mean length of an excursion and Hi the height of excursion i, we have 

Pr (to 4> K)=Pr( sup < 4 ) « Pr f max H, < i] = Pr (Hi < 4)'^/^! . 

\k<K ) \i=\,..\Kli,\ ) 

Thus Pr (ro,4 > K) [\ - Pr (iJi = 4))'^/^l Observing that is of order 1/a', we get that 
the mean value of ro,4 is of order 1/ (aa') = with 

exp I (6.4) Pr (q;q;'to,4 > t) — > e^*/^ as £ gets small. 
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Thus the expected mean time to move from {0} to {4} is the exponential of the global energy 
barrier U (2) — U (0) normalized by e and the time ro,4 normalized by its mean converges 
in distribution to an exponential distribution with mean 2. 

We can check that in the latter model ttq + tti = tts + tt4 showing that the two stable 
state basins share the same weight. 

Suppose the states {0} {2} {4} have a width, say Lq, L2 and L4, where the Markov chain 
undergoes a symmetric random walk before possibly attempting to overcome the energy 
barrier. In this case, the mean values are expected to behave like 

(ro,4) ^ Loi2e(^(2'-^(°»/^ 

(r4,o) ^ L4L2e(^(2)-C/(4))/e 

including a factor involving the characteristic plateaux lengths of the steady states. The 
walker has to overcome its energy barrier but also spends some time in the flat regions 
{0, 2} for the move {0} {4} and {2, 4} for the move {4} {0} . The condition Lq ^ L4 
introduces some skewness in the equilibrium weights of the two stable state basins. These 
considerations are the discrete space-time versions of the result known for a Langevin-type 
stochastic differential equation evolving in a quartic double-well potential U with additive 
white noise with small local variance e. In this context, [IV, if a and b are the stable states 
corresponding to a global minimum of U and if c is the in-between unstable state 

(t h) ~ ^ ^{U(c)-U{a))/e 

^U"{a)\U" (c)| 
^U"ib)\U"{c)\ 



Coming back to the previous symmetric case where {0} {2} {4} are 'simple' states, we 
finally address the following problem: what is the time ro,4 needed to first hit state {4} 
starting from {0} given the walker does not return to {0} again. Note that 

To,4 = T0.4 - sup (fc < To,4 : Xfe = I Xo = 0) . 

We have ro,4 = 1-1- T1.4 where T1.4 is the time needed to first hit state {4} starting from 
{1} of the ergodic chain governed by the transition matrix on {1, .., 4}^: Qi^i ~ 0, Qi^2 = 1; 
Q2.1 = 1/2, Q2.3 = 1/2; (53,2 = a, Q3A = 1 - a; Q4.3 = ct', Q^a = 1 - a'. For this 
Q— chain, the state {1} is now purely reflecting. Using the scale function kp of this new 
chain, Pr = 4) = l/iy9(4) = — > 1/2 (e — >■ 0). Similarly, the mean return time \i to 
state {1} tends to a finite value when e — so that the mean value of to,4 tends itself to a 
finite value when e — > 0. Given there is no possible return to state {0}, the mean time to 
first hit state {4} turns out to be very short compared to T0.4 itself. 
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